Connecting Berry's phase and the pumped charge in a Cooper pair pump 
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The properties of the tunnelling-charging Hamiltonian of a Cooper pair pump are well understood 
in the regime of weak and intermediate Josephson coupling, i.e. when Ej < Ec- It is also known 
that Berry's phase is related to the pumped charge induced by the adiabatical variation of the 
eigenstates. We show explicitly that pumped charge in Cooper pair pump can be understood as 
a partial derivative of Berry's phase with respect to the phase difference cj> across the array. The 
phase fluctuations always present in real experiments can also be taken into account, although only 
approximately. Thus the measurement of the pumped current gives reliable, yet indirect, information 
on Berry's phase. As closing remarks, we give the differential relation between Berry's phase and 
the pumped charge, and state that the mathematical results are valid for any observable expressible 
as a partial derivative of the Hamiltonian. 



Josephson junction devices, e.g. Cooper pair box, 
superconducting single electron transistor (SSET) 
and Cooper pair pump (CPP), have been ex- 
tensively studied both theoretically 1,2 ' 3 ' 4 ' 5,6 and 
experimentallyiiSiSiifliiLiS (For a recent review, see 
Ref . Eii. ) Possible applications include coherent Cooper 
pair pumping 2 with related decoherence studies^ or 
metrological applications^ and the use of these devices 
as superconducting quantum bits (squbits) in quantum 
computationi 1 ! 8 ! 15 In this paper we focus on CPP whose 
idealised tunnelling-cha rging H amiltonian has been 
studied in detail in Refs. I2 ll6ll7l For closed loops in the 
parameter space, we relate the pumped charge to Berry's 
phase, a well-known geometrical phase attained by an 
adiabatically evolving eigenstate of a time-dependent 
Hamiltonianii&iiiSi Some applications of geome trical 
phases in mesoscopic systems, are discussed in Refs. l5l2ll 
and the references therein. We illustrate the results 
both for the SSET and a CPP, and consider the effects 
due to phase fluctuations, present when experimentally 
measuring the pumped current. 

In a CPP the pumping of Cooper pairs is induced by 
cyclic variation of the gate voltages while the evolution 
of the total phase difference across the array, tfi, is ei- 
ther fixed by ideal biasing^ or stochastically decoherent. 
Theoretical predictions are based on the adiabatic evo- 
lution of the eigenstates which splits the induced current 
into two parts The direct supercurrent, which flows 
constantly and is proportional to the ^-derivative of the 
dynamical phase of the eigenstate. The other part, the 
pumped charge, is explicitly induced by the action of 
pumping and proportional to the (^-derivative of Berry's 
phase for closed loops. Existence of such a relation was 
already implicitly stated in Ref. [2J. The underlying rea- 
son for these connections is that the supercurrent opera- 
tor Is is an operator derivative^ of the full Hamiltonian 
with respect to <f>. This also implies that all of the re- 
sults obtained in this paper are valid for any observable 
expressible as a partial derivative of the corresponding 
Hamiltonian. However, in real applications it might be 
reasonable to use the nonadiabatically attained geomet- 
rical phase instead of Berry's phased 



A schematic view of a CPP is shown in Fig. ^ We 
assume that the gate voltages V s j are independent and 
externally operated. The ideally operated bias voltage 
across the array, V, controls the total phase difference, 
(j>, according to d<p/dt = —2eV/h. In the absence of bias 
voltage, <p remains fixed and becomes a good quantum 
numberiSiSi Conversely, the conjugate variable M, the 
average number of tunnelled Cooper pairs, becomes com- 
pletely undetermined. 
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FIG. I: An ideal superconducting array of Josephson junc- 
tions (CPP). Here Ck and Ej t k are the capacitance and the 
Josephson energy of the k th junction, respectively. The total 
phase difference across the array, (f>, is a constant of motion 
in the absence of the external bias voltage V. 



The tunnelling-charging Hamiltonian 
H = H c {q) + H. h 



(1) 



is assumed to be the correct description of the micro- 
scopic system, neglecting quasiparticlc tunnelling as well 
as other degrees of freedom. The charging Hamilto- 
nian Hc(q) depends on the normalised gate charges, 
q := (<7i, . . . , gjv-i); and the number of Cooper pairs 
on each island, n — (m, . . . , njv-i), according to 
(n\Hc(q)\n')(/, = E c (n - q)Sa,H' ■ The function Egjx) 
gives the details of the charging energy, see e.g. Ref. Ho . 
The Josephson (tunnelling) Hamiltonian is given by 



H, = - 



N 

E 

fc=i 



Ej^ k cos(4> k ), 



(2) 



2 



where Ej y k is the Josephson coupling energy of the k th 
junction. The average supercurrent operator can be writ- 
ten in the form 



k=l 



N 



-2e dH 



where the operator derivative^ is defined as 



dH 

7)0 



lim 



H(</>') - H(tf>) 



(3) 



(4) 



In the preferred representation, the parameter space is an 
iV-dimensional manifold IR^ -1 x S 1 , where the elements 
are of the form p = (q, cj>) with (j) G [0, 2Ntt). For N = 
2, the system reduces to a SSET whose Hamiltonian is 
discussed in Ref. 0. Asymptotically exact eigenvectors 
for strong Josephson coupling are given in Ref. |2^. 

We will study the adiabatic evolution of instantaneous 
energy eigenstates \m), while changing the gate voltages 
along a closed path T : t i— > ff(i) = (q (t) , 4>(t)) with 
t £ [0, t]. This induces a charge transfer Q to t '■= Qs + 
Q p , where the pumped charge, Q p , depends only on the 
chosen path, T, determined by the gating sequence, while 
the charge transferred by the direct supercurrent, Q Sl 
also depends on the rate of operation of the gate voltages, 
i.e., the value of r. The total transferred charge, in units 
of — 2e, for state \m) becomes2*i£ 



Qt 



' dc/) 



2 (j) Rc 



(m\M\dm) 



(5) 



Here M = —id/dip is the operator for average number 
of tunnelled Cooper pairs, rj m = — J Q (E m (t)/H)dt is the 
dynamical phase and \dm) is the change in the eigenstate 
|m) due to a differential change dp = (dq, d(f>). A change 
d<p in the phase difference at p induces no pumped charge 
as we find dQ p — 2Im[(M m\M m)]d(f> = 0. In other 
words, the bias voltage V induces no pumped charge for 
fixed gate voltages. 

The expression for Q p is rather similar to the corre- 
sponding Berry's phase 18-19 



7,„(r) = i j> (m\dm). 



(6) 



It should be stressed that Eqs. JSJ) and © are well de- 
fined also for open paths. The derivative d in Eq. © 
is an exterior derivative so, for a closed path T, we may 
integrate Berry's curvature over a two-surface^ 



N N 



T m (r) = iEE 



fe=i i=i 



d(m\ d\m 
Sr dq k dqi 



-dqkAdqi, (7) 



where the boundary of Sp is the path 

We now construct an extended path for which Berry's 
phase is proportional to the charge pumped along the 
path r. Let us define a class of closed paths {T^} by 
L :tH p(t) = (q(t),(f>(t) + ip), where t € [0,r], so 



that To = r. The inverse of a path is the same path 
traversed in the opposite direction, which also holds for 
paths with distinct end points. We define an additional 
class of paths {<pp} according to tpp : t ^ (q,(p + tip), 
where t € [0, 1]. The extended path 



1 cxt 



To ° <pp(p) ° r; 



- 1 o tpZ 1 



p(0) 



(8) 



is also closed and spans a two-dimensional integration 
surface whose width in ^-direction is ip. By traversing 
the boundary the contributions from ^(o) and ^(o) na ^" 
urally cancel, and we find 



T(rr xt )=7(r)-7(r„). 



(9) 



Next we take the limit ^ — > and consider a strip of 
infinitesimal width d(f> between and Tq as illustrated 
in Fig. [3 This means that in Eq. J7J we have either 



FIG. 2: A closed path V has been flattened, i.e. the ends of 
r meet. The infinitesimal pieces of length d<f> are identical, 
but traversed in opposite directions. The two-dimensional 
integration surface Sr,d</> is spanned between orbits V and 
Td4>- Berry's phase corresponding to the boundary of Sr,d<t> is 
identical to — Q v (Y)d4>, with Q P (T) defined in Eq. 1101 . 

dqk = d<j> or dqi = dcj> as the full length of integration and 
we can factor d<f> from the expression for Berry's phase. 
By rephrasing Q p in Eq. 10 as 



N 

Q p (T) = iJ2 f 
k=i Jv 



d(m\ d\m) d{m\ d\m) 
d(j> dq k 



dqu 



dq k , (10) 



we see that it is identical to Berry's phase in Eq. Q 
apart from the factor — d<j>. By taking the limit from the 



equivalent result lirou,. 
first main result 



+0 (7(r) - j(T v ))/tp, we obtain the 



Qs + Q P = -dr) m /d(f) - djm/dcf). 



(11) 



This clearly shows the connection between Berry's phase 
and the pumped charge, which is completely analogous 
with the connection between the dynamical phase and 
the accumulated charge due to direct supercurrent. 

We now proceed in the opposite direction and consider 
strips of finite width instead of infinitesimal ones. By 
integrating the pumped charge with respect to <fi over 
the set {r v }, we obtain the average pumped charge per 
cycle, Q 



p,avc ) 



as 
1 

f Jo 



? p (r )d0 



7(F ) - 7 (rv) 



(12) 



The graphical representation of this situation in a three- 
junction CPP and a SSET, are shown in Fig.|21(I,II) and 
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FIG. 3: I) A projection of path Y in the (51, g2)-plane. II) A 
strip of finite width <p based on the path F. Ill) Fhe fluctua- 
tions of on a single pumping cycle, I. IV) Ideal operation of 
gate voltages produces a strip bounded by the planes 4> — 4>o 
and (j> = 4>o + <p- The same result is obtained approximately 
after many cycles with restricted, stochastic fluctuations of <j>. 

A* 
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FIG. 4: Pumping of Cooper pairs in SSET recuires also an 
ideal control of cf>. The projection of the path T onto g-space 
is a line traversed back and forth. 

in Fig. 4, respectively. The cases are qualitatively differ- 
ent, because there is only one g-coordinate in a SSET. 

We now wish to relate the above results to actual mea- 
surements of Cooper pair pumping. First, consider a 
closed path T corresponding to a fixed value of </>o as in 
Ref. 0. Under ideal operation of gate and bias voltages 
we can change <f> slightly between each cycle and obtain a 
clean strip bounded by the planes <p — 0o and <f> — <fo + V 
as shown in Fig. |3] (IV). Combined with Eq. ijT%|) this 
amounts to an important result for an ideal CPP: The 
measured pumped charge per cycle (i.e. Q p , a vc) yields 
direct information about differences of Berry's phases. 

Obtaining the same information in a real experiment 
is not so straightforward. Neither the phase difference, 
(j), nor the gate voltages are ideally controlled. Never- 
theless, we try to partially circumvent these problems 
using reasonable approximations. First, we assume that 
the gate voltages are operated accurately enough, so that 
the projections onto (/-space nearly coincide. Addition- 
ally, 4> fluctuates stochastically, but these fluctuations are 



restricted during time intervals shorter than the decoher- 
ence time, For times larger than t$ the fluctuations 
mount up too large and the phase coherence of the sys- 
tem is lost. 

The decoherence is induced by any interaction between 
the quantum mechanical system and the environment. In 
a well prepared experiment, e.g. a CPP can be isolated 
from its surroundings so that the main contribution to 
is given by the electromagnetic environment in the vicin- 
ity of the sample and the effects due to finite tempera- 
ture, restricting the measurements to subkelvin regime. 

The decoherence time can be calculated theoretically 
from the fluctuation-dissipation -theorem as in Ref. or 
by looking at coherences, i.e., off-diagonal elements of the 
density matrix. The Hamiltonian in the presence of the 
electromagnetic environment reads 

H = H c (q) + H 3 + H em + H int , (13) 

where H cnv = Y^jiK^j + l/2)7iWj and bj and bj are the 
creation and annihilation operators of the bosonic envi- 
ronmental mode j with energy hujj , respectively^ As an 
example, we consider a SSET but it should be stressed 
that the result generalises for any number of junctions. 
We write the density matrix p k = ^> fc vfi fe t m the basis 
consisting of two SSET charge states, {|m)}^ =0 and en- 
vironmental modes {\k = (k\, fc 2 , . . .))}. Then the Hamil- 
tonian describing the interaction between SSET and the 
environment becomes2Z*2& 

a.* = -iV*E ^VJ ~ 6 9 ( -1 ) ' (14) 

where Zj is the impedance of the mode j and Rk = 
h/e 2 ~ 25.8 kfi is the resistance quantum. 

The equation of motion for p\ in the interaction pic- 
ture is given by the Liouville equation, iK(dpj (t) / dt) = 

[-Hint, 1 ■ Pi(t)}. By solving the differential equation for the 
coherence matrix elements and tracing out the environ- 
mental configurations k we obtain the final result 

01,12(0 = Pi,i2(0) exp [-2ReJ(i)] , (15) 

which corresponds to the same time scale as given 
by the fluctuation-dissipation -6-theorem£ Here J(t) 
is the phase-phase correlation function J(i) = 

([ip(t) — 95(0)] ip(0)} and pi = (pi)% -~ I n case of purely 
resistive electromagnetic environment, R e , Eq. i|15|) 
yields^ ~ [7r^/(16fcBT)]i?K/-R e , where we have as- 
sumed nonzero temperature and irk^Tt/h ^> 1. For re- 
alistic measurement parameters, e.g. T = 10 mK and 
R c = 10 f2, one obtains a rather long time r$ ~ 0.4 ps. 

Returning to Berry's phase, we assume an initial value 
4>o and consider time intervals shorter than , effectively 
restricting <j) to a finite range [4>i, 4>2\ 9 4>o- If sufficiently 
many (identical) cycles of gate voltages are performed 
during this time, the fluctuations of 4> yield a relatively 
thick mesh of trajectories within the strip. Although 
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the 'weights' for different values of (j> are uneven, we ap- 
proximate the mesh with a uniform distribution which 
is a subset of the range [4>i,4>2]- This corresponds to a 
well-defined strip as in the ideal case of Eq. 1)1 2|l and is 
presented in Fig. [3] (IV). A cycle, £, with exaggerated 
fluctuations in cf), is shown in Fig. |21 (HI). Due to the 
stochastic nature of the fluctuations, it is impossible to 
predict the correct range to be used. Nevertheless, for 
periods that are short enough, the correspondence be- 
tween Berry's phase and the measured pumped charge 
exists in the sense of Eq. 1)12)1 . If the end points of the 
full pumping cycle are sufficiently close, Eq. is valid, 
at least in the framework of the model. 




»d<f> => 



pumping 

dq 
« i 



Ym(Q 



equivalent Berry cycle 



consider Berry's phase j m induced by an infinitesimal 
closed cycle C at (q, <p) with sides dq and d(j) as shown by 
the left-hand-side of Fig. On the right-hand-side, the 
pumped charge due to dq multiplied by d(f>, is identical to 
Berry's phase induced by the discontinuous path below it. 
By following any closed pumping path T and integrating 
the pumped charge, we recover Eq. If the path is 

not a closed one, a nontrivial integration with respect to 
<fi remains, regardless of the width of the strip. 

In conclusion, we have shown explicitly how the 
pumped charge in Cooper pair pump can be understood 
as a partial derivative of Berry's phase with respect to 
the phase difference cf> across the array. We have only 
used the fact that the supercurrent operator Ig is an 
operator derivative of the full Hamiltonian. Thus these 
results generalise for any observable with this property. 
We have also shown how one could obtain information 
about Berry's phase by measuring the pumped current 
in a CPP. 



FIG. 5: An infinitesimal cycle C corresponding to Berry's 
phase 7m (C) consists of four legs. The charge transfer Q p for 
a fixed (j) is identical to Berry's phase induced by traversing 
the legs in the directions shown on the right-hand-side divided 
by d<j>. This path can not be continuously followed in the 
parameter space. 

As a final note, we construct the differential relation 
between Berry's phase and the pumped charge. Let us 
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